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abstract 

We propose a mechanism based on the rubble-pile hypothesis of the cometary nucleus (Weiss- 
man 1986) to explain the catastrophic breakup of comet LINEAR (C/1999 S4) observed during 
July- August 2000. We suggest that a solid nucleus made up of 10-100 m “cometesimals” (Weiden- 
schilling 1997) contains a network of inter-connected voids in the inter-cometesimal regions. The 
production of super-volatile (i.e., species more volatile than water) gasses into these voids occurs 
due to the thermal wave propagating through the nucleus and associated phase transitions of wa- 
ter ice. The network of voids provides an efficient pathway for rapid propagation of these gasses 
within the nucleus resulting in gas pressure caused stresses over a wide regime of the nucleus. This 
provides a mechanism for catastrophic breakups of small cometary nuclei such as comet LINEAR 
(C/1999 S4) as well as for some observed cometary outbursts including those that occur at large 
heliocentric distances (e.g M West et al 1991). We emphasize the importance of techniques such 
as radar reflection tomography and radiowave transmission tomography (e.g., Kofman et al 1998) 
aboard cometary missions to determine the three dimensional structure of the nucleus in particular 
the extent of large scale voids. 
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1. Introduction 

Comet LINEAR (C/1999 S4) was originally discovered by The Lincoln Near Earth Asteroid 
Research (LINEAR) program of the MIT Lincoln Lab on September 27, 1999 as an asteroidal object 
whereas the cometary nature of the object was noted by Durig et al (1999) a few days later. The 
close perihelion distance (0.77 AU) of the orbit and the prediction of a likely naked-eye visibility in 
July 2000 (Green 1999a, 1999b) attracted attention to this comet. As a result, it was monitored by 
different observers during the following months leading to the perihelion passage. The observations 
around the time of perihelion passage provided the details of the unexpected catastrophic breakup 
of comet LINEAR (C/1999 S4). These were the first such observations of a complete disruption of 
a cometary nucleus while the disruption was taking place (Kidger 2000, Licandro 2000, Weaver et 
al, 2001), providing an insight into the death of a comet, and in the process revealing clues about 
its structure. 

In this paper, we propose a model to explain the catastrophic breakup of comet LINEAR 
(C/1999 S4). The model will additionally provide a viable mechanism which can explain some of 
the cometary outbursts for other comets — especially outbursts originating at large heliocentric 
distances (and also from the nightside). 


2. Comet LINEAR (C/1999 S4) 

Marsden et al. (1973) introduced radial and transverse nongravitational parameters, A\ and 
A 2 , with the nongravitational acceleration of the comet, a ng , due to outgassing forces expressed by 
(neglecting the component normal to the orbital plane) 

a ng = ^1 g{rh)i + A 2 g{r h )t (1) 


l The National Optical Astronomy Observatory is operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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where T} x is the heliocentric distance and g{r^) describes the outgassing rate as a function of the 
heliocentric distance, f is the instantaneous orbital unit vector in the outward radial direction while 
t is the instantaneous transverse unit vector directed normal to r, in the orbital plane, and in the 
direction of the orbital motion. 

Astrometry of comet LINEAR (C/1999 S4) revealed that it experienced comparatively large 
nongravitation al accelerations when compared with those of other comets (cf. Marsden and Williams 
1993); for LINEAR (C/1999 S4) Ai=+9.5±0.23 and A 2 =-2.11 ± 0.09 in units of 10“ 8 AU day' 1 
(Marsden 2000). Even for the highly active comet IP/Halley with localized outgassing (Keller 
et al 1986), the corresponding nongravitational parameters are smaller by more than an order of 
magnitude. As nongravitational accelerations are directly proportional to the net outgassing forces 
and inversely proportional to the nuclear mass (e.g., Yeomans and Chodas 1989, Samarasinha and 
Belton 1995), these large values of nongravitational parameters suggest a comparatively small nu- 
cleus. Our prior experience with modeling nongravitational forces (e.g., Samarasinha and Belton 
1995, Samarasinha et al. 1996) indicates that the nuclear radius is of the order of one kilometer. 
We note that the exact number depends critically on the details of outgassing. 

An independent lower limit to the nucleus radius can be derived from the total water production 
rate. Schleicher and Eberhardy (2000) derive a water production rate of 2.6xl0 28 molecules s" 1 
(based on observations of OH) on July 13, 2000 when the heliocentric distance to the comet was 
0.81 AU. Simple thermal models (e.g., Cowan and A’Hearn 1979, Julian et al. 2000) indicate 
this corresponds to a minimum nuclear radius of % 0.7 km provided water was produced only by 
nucleus outgassing from the sunward hemisphere and the nucleus was not in an outburst during 
the observations. 

Above arguments, while not conclusive, strongly favor that the nuclear radius of comet LIN- 
EAR (C/1999 S4) is of the order of one kilometer. Therefore, it is smaller than the recent bright 
comet, Hale-Bopp (C/1995 Ol) (radius of 20-35 km; Weaver and Lamy 1997), and another exten- 
sively observed comet, IP/Halley (effective radius of about 5 km; Keller et al. 1986). 

Weaver et al (2001) estimate that the largest breakup fragments of LINEAR (C/1999 S4) are 
consistent with diameters of about 100 m. This is consistent with the suggestion by Weidenschilling 
(1997) that 10-100 m size cometesimals are the likely building blocks of cometary nuclei based on 
his numerical simulations of growth of icy planetesimals in the outer solar system. 


3. Splitting/Breakup of Cometary Nuclei 

The splitting events, primarily the breaking away of cometary fragments from the nucleus 
are not uncommon. Based on the observations of 49 comets, Chen and Jewitt (1994) derive a 
lower limit of 1 comet splitting per comet per 100 years. There is extensive literature on the 
split comets including a review by Sekanina (1982). In a more recent paper, Sekanina (1997) 
revisits the subject and prefers a characterization based on whether the splitting occurred due 
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to tidal interactions or not. Observational properties of split comets seem to depend on whether 
the splitting mechanism is of tidal origin or not. In particular, tidal disruption tends to produce 
fragments in a line, whereas other modes of disruption produce scattered fragments, as was the 
case for comet LINEAR (C/1999 S4). 

There are a number of breakup mechanisms discussed in the literature for non-tidal splittings. 
These include rotational splitting (e.g., Sekanina 1982), including rotational splitting facilitated by 
evolution of the nucleus into a more elongated shape due to volatile loss (Hartmann and Tholen 
1990), rotational splitting and/or gas pressure caused breakup facilitated by cracks of the nucleus 
which may occur either during phase transitions (e.g., Smoluchowski 1981, Jenniskens and Blake 
1996) or via thermal stresses (e.g., Iviihrt 1984, Tambovtseva and Shestakova 1999). 

Some of the mechanisms proposed to explain cometary outbursts can also lead to cometary 
breakups like the impact of an asteroidal object (e.g., Hughes 1991, Kadano 1999) and the exother- 
mic polymerization of HCN (e.g., Rettig et al 1992). 

A model for explaining the breakup of the nucleus of comet LINEAR (C/1999 S4), unlike 
that for isolated splitting events, should explain why the whole nucleus underwent a catastrophic 
breakup. In other words, the model should provide a mechanism which can break the whole nucleus 
into multiple fragments abruptly consistent with the observations. 

Rotational splitting is most likely to occur near the “equator” of the nucleus where the dis- 
ruptive force due to rotation is the largest. It does not explain how the entire nucleus can be 
broken apart. A collision with an asteroidal object is most likely to occur near the ecliptic plane. 
However, the breakup of comet LINEAR (C/1999 S4) occurred close to perihelion when the comet 
was well away from the ecliptic plane. Furthermore, numerical studies by Asphaug et al (1998) 
indicate energetic collision of a house size asteroidal object with a kilometer-sized rubble pile target 
is unlikely to disrupt the target. Finally, existing gas pressure based breakup mechanisms do not 
explain how the entire nucleus can be involved in a catastrophic breakup. 


4. Model 

We propose that the solid nucleus is made up of an aggregate of 10-100 m “cometesimals” 
(Weissman 1986, Weidenschilling 1997), and contains significant voids in the inter-cometesimal 
regions (see Fig 1 for a cartoon of the nucleus model). This network of inter-connected voids 
provides an efficient pathway for rapid propagation of super-volatile (i.e., species more volatile 
than water — dominated by CO and C0 2 ) gasses over a wide regime of the nucleus causing gas 
pressure related stresses. These stresses, if sufficiently large, can eventually breakup the nucleus 
into the original “cometesimals” (i.e., its original building blocks). 

The super-volatile gas production into inter-cometesimal voids can occur due to the inward 
propagation of the thermal wave as the comet comes closer to the sun. The exothermic phase 
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transition of water ice from amorphous to crystalline may further facilitate the efficient inward 
propagation of the heat front (Prialnik and Bar-Nun 1987). Super- volatile gas release is associated 
with this phase transition (Prialnik and Bar-Nun 1987), due to restructuring of the water ice 
matrix caused by transitions between different forms of amorphous water ice occurring at lower 
temperatures (Jenniskens and Blake 1994), and possibly due to the exothermic polymerization of 
HCN (Rettig et al 1992). 

The propagation depth (skin depth), 6 , of the heat front is given by 


S = 



I t ( T \ l t 2 
~pC \n) 


( 2 ) 


where AT, p, C, and It are thermal conductivity, density, specific heat, and thermal inertia of 
cometary ice respectively, while r is the period of the heat wave. For a cometary-ice thermal 
inertia of 0.016 cal cm -2 s _ 2 K -1 (which is based on model fits to water production rates in comet 
lP/Halley; Julian et al. 2000), over the pre-perihelion leg, the propagation depth of the heat front 
for comet LINEAR (C/1999 S4) is of the order of 10 m. The exothermic phase transitions described 
earlier would further enhance the propagation depth of the heat front. Comparable depths were 
derived by Prialnik and Bar-Nun (1987) for the crystallization front in case of a comet in the orbit 
of IP/Halley. Likewise, Benkhoff and Huebner (1995) using an entirely different thermal model for 
the depth of the heat front for a comet in the orbit of 73P/Schwassmann-Wachmann 3 obtained 
similar results. Therefore, we argue that the thermal wave had sufficient time to propagate one 
cometesimal width (at least 10 m but likely more), in the case of comet LINEAR (C/1999 S4) by 
the time it broke up near perihelion. 


The gravitational potential energy, E g , of a constant density spherical nucleus with radius Rn 
is given by 

3 GM 2 l ^ 2 - G p 2 R N 5 ( 3 ) 


Eg 5 R n 


15 


where M is the mass and p is the bulk density of the nucleus while G is the universal gravitational 
constant. Therefore, for a nucleus of Rn ^1 km, E g «10 18 ergs. Now, by assuming zero tensile 
strength against breakup, one can derive the minimum super-volatile gas mass required to com- 
pletely disperse the nucleus due to gas pressure. The average kinetic energy for a mole of gas, Ek , 
is given by 


E k = 5 RT 


( 4 ) 


where T (p^IOO K) is the temperature and R is the universal gas constant. From the above two 
equations, to completely disperse the nucleus due to super- volatile gas pressure, we derive a min- 
imum gas mass of 10 8 moles. This gas mass (when in the form of ice) would occupy only a small 
fraction (of the order of 10 -6 ) of the volume of a Rn ~1 km nucleus and is consistent with our 
model. The solar constant (i.e., incident solar energy at 1 AU) is 1.4xl0 6 ergs s'" 1 cm -2 . Therefore, 
even if one assumes that the ultimate energy source for the breakup of the nucleus is solar energy 
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(actually in the model, it acts only as the catalyst), a very small fraction of the total incident solar 
energy over the pre-perihelion leg can provide the necessary breakup energy. 

Thermal models of the nucleus (e.g., Benkhoff and Huebner 1995, Capria et aL 1996, Enzian et 
al 1999) which adopt a porous media for cometary ice find efficient gas diffusion through a series of 
unhindered tubal pores of radii typically of the order of mm or less. In those models where super- 
volatiles get depleted in the top most layers, they are produced at a depth and transported through 
pores in the porous media. Therefore, how can gas pressure buildup occur in the inter-cometesimal 
voids in our model? We suggest, that the process of collisional accretion leading to the formation 
of 10-100 m cometesimals (Weidenschilling 1997) did alter the structure of the accreting taget 
bodies. This structure alteration includes compaction and local melting at the collision sites — a 
sentiment expressed in Weissman (1986). Collision caused compaction of porous material is argued 
in Sirono and Greenberg (2000) and also seen in collision studies involving porous targets aimed 
at explaining large craters of the asteroid Mathijde (Housen et al. 1999). It is highly unlikely that 
cometesimals formed via a series of collisions could retain a set of unobstructed pores extending up 
to tens or hundreds of meters depth. When the porous media is modeled by a series of tubal pores, 
introduction of an effective scale length, /, for the “unobstructed pore length” is needed to account 
for the structure alterations discussed above. Unfortunately, to this author’s knowledge, there are 
neither numerical nor laboratory studies on structure alterations applicable to icy cometesimals 
to derive a reliable number. While cometesimals can be porous, it is unlikely that their structure 
allows unhindered diffusion of gasses over tens or hundreds of meters. For example, for a depth of 
10 m, l ^1 m would result in a decrease of more than four orders of magnitude (i.e., e -10 ) for the 
effective number of available “unhindered” tubal pores. This facilitates gas pressure build up at 
high depths and has an effect similar to what a mantle does in the model of Mohlmann (1996). 

Further evidence for our current insufficient understanding of the porous structure of cometary 
material and gas flow through it is apparent from the inability of the existing thermal models to 
accurately predict the observed CO production rate of comet Hale-Bopp (C/1999 Ol) as a function 
of the heliocentric distance (Bockelee-Morvan and Rickman 1997). Therefore, further refinements 
to current thermal models, especially gas transport through cometary material, as well as exploring 
new mechanisms of super-volatile transport through the nucleus are required. 

The other condition which prohibits a pressure build up in the inter-cometesimal voids is the 
presence of unobstructed large void outlets to the coma. Clearly, LINEAR (C/1999 S4) which broke 
up catastrophically (which was a rare event) was different from the majority of comets which do not 
undergo similar demises. We suggest what differentiates LINEAR (C/1999 S4) from other comets 
is lack of such void outlets to the coma thereby facilitating internal pressure build up. Blockage of 
surface outlets can occur due to reaccretion of “regolith” during the collisional accretion process 
(for example, see Leinhardt et al 2000) and/or due to regolith redistribution associated with 
sublimation induced surface erosion — a surface process which has been proposed earlier (e.g., 
Sekanina 1991). For the purpose of calculating the rate of super-volatile gas effusion from any 
residual narrow openings to the coma at the blocked void outlets, let us define A to be the sum of 
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cross sections of all such openings. From kinetic theory, the relaxation time, r outl of the confined 
gas is given by 

V f2wfj, 

Tout= A\~RT ^ 

where V is the volume of the gas (assumed to be equal to the volume occupied by inter-cometesimal 
voids) and f.i is the molar weight of gas. 


In our model, the super- volatile gasses produced into voids can quickly spread throughout the 
network of inter-cometesimal voids. From the application of Bernoulli’s principle and the ideal gas 
law, we derive 


v 


I2RT 




( 6 ) 


where v is the gas propagation speed in the voids. Therefore, the propagation of gasses in inter- 
cometesimal voids occur on a timescale, Ti nt) given by 


2 Rn 

Tint & 

V 


(V 


(i.e., of the order of 10-100 s). For the confinement of super-volatiles and the subsequent pressure 
buildup, T oui > Tint * he., 



(8) 


where /(= ^ ^ -3 <C 1) is the fractional volume of the inter-cometesimal voids. In the case of 

£ It fl-PJ 

LINEAR (C/1999 S4), A must be < 10 9 cm 2 (i.e., less than 1% of the total surface area). 


For the inter-cometesimal voids closest to the surface, the temperature of the cometesimal 
boundaries can be approximated by the amorphous to crystalline water ice transition temperature 
(^ 137 K), but the cometesimal boundaries further inside the comet are likely to be much cooler, 
unless there is an interior energy source such as radiogenic heating under suitable conditions (e.g., 
Prialnik and Podolak 1995 and references therein). Consequently, condensation of super- volatile 
gasses can occur if P x > P/(T), where Pi is the gas pressure and Pf(T ) is the saturation gas pressure 
for a super- volatile species i. Therefore, the upper limit to the total gas pressure, P uvvtr , can be 
expressed as 

Puppe r ='£P°(T). (9) 

i 

Breakup of the nucleus would occur if the disruptive forces (i.e., forces due to gas pressure and 
rotation) exceeds the binding forces due to self-gravity and large-scale tensile strength which holds 
the cometesimals together. We can mathematically express this as (see Appendix) 


_ 27 t 2 oRn 2 cos 2 A ^ 2 ^ n _ o 

P + fLJL ^ > -7T Gp 2 R N 2 + a 

TN 2 3 


( 10 ) 


where P is the total gas pressure, rjv is the rotation period of the nucleus, A is the latitude, and a 
is the large-scale tensile strength. By rearranging this equation and by substituting P by P upp e r, 
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we can derive an upper limit to the large-scale tensile strength. 
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For illustration purposes, by assuming that the gas pressure is dominated by the saturated gas 
pressure P^ Q of the most volatile and abundant super-volatile CO and neglecting the term due to 
rotation (even on the equator, rotation term can exceed the self-gravity term only if the rotation 
period is smaller than ^2 hrs where p is expressed in terms of g cm" 3 ; also Pravec and Harris 
2000), we have 

G upper — P S * * * * CO “ g KGp 2 RN 2 . ( 12 ) 

Saturated gas pressure of CO can be approximated by 

PZo = Acoex p(-^p) ( 13 ) 


where Acro = 1*2631 x 10 10 * * dyne cm -2 and B(%)=76 4.16 K (Fanale and Salvail 1990). Fig 2 depicts 
the dependence of a upper as a function of P£o and R n. The binding effect due to self-gravity 
becomes relevant for a Rn =1 km nucleus if a < 10 3 dyne cm -2 while for a Pjv =10 km nucleus this 
holds true if a <10 5 dyne cm -2 (i.e., about 10 -3 and 10 -1 atmospheres respectively). Therefore, 
self-gravity can play an increasingly larger role in preserving larger nuclei. 


5. Discussion 

A major feature of the proposed model is the ability of super-volatile gasses to rapidly prop- 
agate within the inter-cometesimal voids. These natural voids created during the formation of 
a cometary nucleus in the solar nebula provide an efficient pathway for gas propagation. In the 
rubble pile model for cometary nuclei (Weissman 1986), the tensile strengths at which constituent 

“cometesimals” held together are lower than the tensile strengths of individual “cometesimals” 
themselves, thus providing a set of natural boundaries for breakup. Such a scenario is consistent 
with observations where the nucleus of comet LINEAR (C/1999 S4) broke into fragments of whose 
sizes are approximately in the same range as those predicted by Weidenschilling et al. (1997) for 
the basic building blocks of comets. 

The proposed scenario raises the issue of why subsurface gas pressure broke up comet LIN- 
EAR (C/1999 S4) whereas not other comets. Inter-cometesimal voids lacking sufficiently large 

unobstructed outlets to the coma which facilitates internal gas pressure build up, as we suggest for 

LINEAR (C/1999 S4), differentiate it from the majority of comets which do not undergo similar 
catastrophic breakups. In addition, small cometary nuclei such as comet LINEAR (C/1999 S4) are 

more prone to such complete catastrophic breakups due to gas pressure driven mechanisms than 

larger nuclei. Larger cometary nuclei are more stable because of the combined effects of (a) the 

ability of the larger self-gravity to hold the nucleus together (which scales as Rn 2 ', see previous 
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section), (b) a larger ratio of volume of voids to rate of super- volatile production (which scales 
as Rjy ), and (c) the increased likelyhood of voids having a large unobstructed outlet to the coma 
(which scales as Rn 2 )- 

In a different context, large scale subsurface propagation of gasses was previously proposed by 
Mohlmann (1996), where he suggests that the active regions of water outgassing on cometary nuclei 
are simply local “source sites” of accumulated subsurface flows from distant regions in the nucleus, 
in contrast to water sublimation occurring predominantly at the surface in the conventional thermal 
models. Model results from a scenario as proposed by Mohlmann (1996) need to be compared with 
observations, in particular the water production rate and its variations with heliocentric distance 
and nucleus rotation. 

If the inter-cometesimal voids have an unobstructed outlet to the coma, this same physical 
mechanism can be responsible for some of the observed cometary outbursts, such as those observed 
at large heliocentric distances. A likely example is the 14 AU outburst of comet 1 P/Halley (West et 
ai 1991). Observationally verifying whether an outburst is related to a super- volatile gas pressure 
driven mechanism is possible as there should be enhanced abundances of species more volatile than 
water during an outburst. In fact, such an association has been made by Feldman et ai (1986) in 
the case of comet lP/Halley. 

The issues raised in this paper point to the necessity of detailed numerical and laboratory 
studies of porous structure of collisionally formed “cometesimals” as well as detailed thermal models 
of the nucleus based on inputs from such studies in the future. 

Finally, we wish to point out the importance of effective three dimensional structure deter- 
ministic techniques such as radar reflection tomography (M. Belton, private communication) and 
radiowave transmission tomography (e.g., CONSERT experiment aboard Rosetta; see Kofman et 
al. 1998) in future cometary missions. Such determinations would not only tell us the three di- 
mensional structure of cometary nuclei, but also be useful to assess the likely role of subsurface 
structure in cometary activity and breakup events. 

Appendix 

The inequality for the non-tidal breakup of a nucleus can be expressed as P + Fr > Fq + <? 
where P is the total interior gas pressure, Fr is the disruptive force due to rotation per unit area, 
Fq is the binding force due to self-gravity per unit area, and a is the large-scale tensile strength 
holding cometesimals together. We derive Fr and Fq for a spherical nucleus. 

Derivation of Fr: 

Consider a ring of cometary material at latitude A having width rdX and height dr. The 
disruptive force due to rotation is given by mrtr ^ - v ~ = ( 27rrp cos V T ^) where m r i ng is 
the mass of the ring, V is the speed of rotation, r is the distance from the center of the nucleus, 
p is the bulk density of the nucleus, and is the rotation period. As this force acts upon an 
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area of 27rrcosA rdA, the force per unit area is given by 4rr p - c °^ Xr dr . By integrating we find 

Ffi = 27T - across the entire nuclear radius R n. This expression for Fr is consistent with 

the one cited in Sekanina (1982) for A = 0. 

Derivation of Fq: 

Consider a spherical shell of cometary material having thickness dr . The gravitational force on 
the shell is given by Gm ^^sh.Qii — )^ 47 l, pr ,^ r ). where G is the universal gravitational constant, 

rrii n t is mass interior to the shell, and m s heii is mass of the shell. Since this force is acted upon an 
area of 47rr 2 , the force per unit area is given by |t xGp 2 rdr. By integrating we get Fq = | nGp 2 R n 2 
across the entire nuclear radius. 
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Fig. 1. — A cartoon showing a 2-D cross section of the 3-D nucleus model containing inter- 
cometesimal voids and individual cometesimals. In the cartoon, the cometesimals are shown in 
light gray while arrows denote regions where outgassing can occur. Outgassing of volatiles (domi- 
nated by water) can occur from non-mantled cometesimal surfaces. In addition, super-volatiles can 
escape from void outlets not covered by regolith/debris. In the case of comet Linear (C/1999 S4), 
we suggest that the regolith/debris obstructs all the void outlets, thereby facilitating internal pres- 
sure buildup. Note that the cartoon is not drawn to scale (for example, the number of cometesimals 
is underrepresented and small scale features are exaggerated). Therefore, we caution the reader 
not to overinterpret the cartoon. 

Fig. 2. — Dependance of <J upper as a function of P^ 0 (as well as T) and R^j. Nuclear bulk density 
of 0.5 g cm -3 is assumed in calculations. For a given nuclear radius, unless the tensile strength is 
sufficiently large, the disruptive forces can break/split the nucleus as self-gravity cannot hold the 
nucleus together. This critical tensile strength increases with nuclear radius. 
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